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The transistor as an active two-port network 


by A. R. Boothroyd*) 


Compiled from invited lectures on research on transistors given at the Institute for Advanced Electrical 
Engineering at the Swiss Federal Institute of Technology, Zurich, on July 5—7, 1960 


Der Transistor als aktiver Vierpol. In dieser Arbeit werden fundamentale Beziehungen für den Transistor- 
vierpol betrachtet. Zuerst werden die Bedingungen für dynamische Stabilität bei beliebigen Quellen- und 
Lastimpedanzen mit Hilfe des Sfernschen Stabiitätsfaktors k behandelt. Es wird gezeigt, dass k weder 
einen allgemeinen Ausdruck für einen Stabilitäts- oder Instabilitäts-«Spielraum» gibt, noch ein prin- 
zipielles Mass des Übertragungsgewinnes bei Stabilität darstellt. Der Grund dafür ist, dass für gegebene 
Quell- und Lastimpedanzen k je nach den benützten Netzwerkparametern (h, z, y oder g) verschiedene 
Werte annimmt. Ës wird nun dargelegt, dass ein Faktor n für die Beschreibung des Übertragungsge- 
winns geeignet ist, umsomehr als letzterer für starke Fehlanpassung am Ein- und Ausgang des Transistors 
eine vom gewählten Parametersystem nahezu unabhängige Grôsse wird. Mittels des Faktors # wird eben- 
falls der optimisierte maximal verfügbare Leistungsgewinn (MAG) eines unilateralisierten Netzwerks aus- 
gedrückt. Im darauffolgenden Kapitel wird ein vom Parametersystem vollständig unabhängiger Faktor S 
definiert, welcher unter bestimmten Umständen in k bezw. n übergeht; dieser Stabilitätsfaktor S ist 
deshalb eine allgemeinere Grôsse als die Faktoren k oder #. In einem weiteren Abschnitt wird die Berech- 
nung von Rückkopplungsnetzwerken für «optimale» Instabilität von Transistor-Oszillatoren behandelt. 
Auslegung auf maximalen negativen Ausgangsleitwert ergibt dabei Schwingungen bis zu Frequenzen, wo 
der MAG des Transistors gleich Eins wird, oder für maximalen Belastungsleitwert bei beliebiger Frequenz. 
Ein beträchtlich verändertes Rückkopplungsnetzwerk resultiert dabei aus der Forderung nach Verfügbar- 
keit der maximalen Klasse À — Ausgangsleistung. Schliessiich werden noch die Bedingungen für eine 
gute Frequenzstabilität von Transistor-Oszillatoren untersucht. 


Properties of the transistor as an active two-port network are presented in fundamental terms. First, the 
conditions for dynamic stability, or otherwise, with arbitrary terminations at the ports, are considered in 
terms of Sfern’s stability factor k. It is shown that k does not give a basic expression for à “margin” of 
stability or instability, nor is this parameter a basic measure of transducer gain performance when the 
network is stable. The reason for these properties is that for given terminations k takes different values for 
different sets of network parameters (h, z, y or g) that might be used for analysis purposes. À performance 
factor n is shown to be appropriate for the expression of gain performance, especially since for a condition 
of large mismatch at the ports of the transistor the transducer gain then closely approaches an invariant 
with change of parameter system. The optimised MAG of unilateralised amplifier networks is also ex- 
pressed in terms of #. An exact invariant factor S (i. e. independent of the parameter system) is defined, 
which under special conditions is seen to take on the roles of k or #; this invariant factor is thought to be 
more basic than the other two parameters. The design of fecdback networks to cause ‘‘optimised” 
instability in transistor oscillators is discussed. Optimisation for maximum negative output conductance 
at a port ensures oscillation at frequencies up to where the MAG of the transistor falls to unity, or for 
maximum load conductance at any frequency. À considerably different feedback network is shown to 
result if, alternatively, it is required that the maximum class-A power output be developed. Finally 
considerations relevant to the realisation of good stability of oscillation frequency are examined. 


*) This is the first of two papers by Dr. À. R. Boothroyd on recent work on transistor theory carried 
out by his research group at Imperial College of Science and Technology, London. 


1. Introduction 


The transistor is an active three-terminal circuit element; it is in general bilateral in any 
arrangement and is potentially unstable over a wide frequency range. In each of the three useful 
two-port amplifier configurations —common-base, common-emitter and common-collector— 
the maximum available power gain (MAG), achieved with conjugate matched terminations, is 
finite over restricted low and high frequency ranges, as illustrated in fig. 1. Between the frequen- 
cies fa and fp, the values of which depend upon the transistor configuration, the concept of MAG 
is meaningless; with suitable terminations any transducer gain up to infinity may be realised, 
or the circuit may be dynamically unstable. The limiting frequency fmax at which MAG is 
unity may be shown to be invariant with change of device configuration [1]. This is the maximum 
frequency of usefulness of the transistor as an amplifier, and the maximum frequency of oscil- 
lation with a lossless external feedback network. 


0 
MAG Û 
| range of 
dB (2 ntial -,1 
Û instability | 
Û 
i Û 
! 
Fig. 1. 
© Typical dependence of maximum available gain of a transistor 
£ f Re on frequency. 


For use of the transistor as an amplifier over the frequency range 0 f< fmax stability may 
always be achieved by either 


(a) suitably mismatching at the ports 


(b) applying external feedback. 


These measures may be applied separately or in combination such that an adequately large well- 
defined power gain is realised. If in (b) the reverse transfer parameter of the transistor is reduced 
to zero by applied feedback, the overall two-port network is then unilateral. In the design of 
transistor amplifiers a number of questions arise: viz. which device configuration is most 
suitable; by which means should dynamic stability be attained, in the case of a potentially unstable 
network; if unilateralising feedback is to be adopted, what circuit arrangement scheme should be 
used for best gain performance ? The main purpose of the present paper is to indicate how these 
and related questions may be answered analytically. 

Schemes of external feedback may be classified in terms of the corresponding systems of 
two-port parameters most appropriate for their analysis, or «parameter environments». The 
performance of à two-port network may be analysed in terms of any of six alternative sets of 
parameters; of these possibilities four sets—h, z, y and g—are most commonly used, and they 
define the four basic feedback schemes indicated in fig. 2. In this figure, the source and load 
terminations of the overall two-port network are series- or shunt-expressed (1. e. in impedance or 
admittance form) according as current or voltage respectively is the independent variable for the 
chosen parameter system at the port concerned. Given certain internal parameters for the active 
two-port network (transistor) each of the passive feedback systems may be chosen to render the 
overall (combined) network unilateral; the feedback may further be optimised so that the greatest 
MAG is realised. In general the values of optimised MAG attained in the four feedback «environ- 
ments» are different, and one particular system gives the best performance. 
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Fig. 2. The four basic feedback «environments». 


9 - Environment 


External feedback may also be applied to the transistor in order to cause dynamic instability, 
as required in the case of an oscillator. Again, in the design of transistor oscillators a number of 
questions arise : viz. how should the feedback system be arranged to secure the greatest instability— 
greatest regeneration—at à given frequency; under what conditions does an oscillator deliver the 
maximum power output; what are the requirements for good frequency stability ? Answers to 
these questions are given in the following sections. 


2. Stability and gain performance 


2.1. Sfability 
Analysis of the conditions for the dynamic stability, or otherwise, and of the gain performance 
of an active two-port network, may proceed in terms of any of the possible sets of two-port 


parameters (say h, z, y or g). In general, let the selected parameter set be designated p as in 
fig. 3 where 


P—p + jo —h = yon, 


and the source and load terminations ps, pp are written in the form 2 = r + jx or y — g + jb 
as appropriate. Further, let the total self parameters at the two ports be written 


p1 = 01 +]J01—=Ps + pui 


pa = pa +] 62 — pr + pos 


Fig. 3. 
Definition of generalised 2-port parameters and terminations: 
p = h,2,y, Or g. 


Total self parameters : P, = b.+b, 
mem | 

Total port paramekers  : P, = P+P. 
Pre RP 

In general : P = p+jr 


Then, according to Stern [2] the network is absolutely stable 1f the stability factor 


(1) where pa1 p12 = M + jN 
and | pa pr2| = L = PM? + N? 


The term «absolutely stable» is in the sense that for given real parts of terminations ps, px, the 
network is stable for any 6, 67,; a condition that is particularly relevant in the case of an amplifier 
which is to be adjusted for «tuning» at the ports. Thus stability may be achieved either by 
adequate mismatching at the ports (to make p1, 02 sufficiently large), or by applying feedback 
such that p19 is suitably small, and possibly reduced to zero to achieve unilateralization. 

An important property of the stability factor k is that in general for k 1, 


kn£k:£ky£ ke 


for given physical terminations (suitably expressed according to the parameter system p); thus 
k does not express an absolute or invariant «margin» of stability or instability. Moreover, a given 
value of k > 1, 1.e. 

ER ERA 


defines different terminations and values of transducer power gain in the different parameter 
environments, so that k is not an absolute or invariant expression of gain performance. Only 
under the special condition k — 1, defining the boundary between stability and (possible) 
instability, is k an invariant. 

A useful network parameter, that characterises the active two-port network, is the «inherent 
stability factor» 


(2) R— 2P11022 
L+M 


this is the value of k for zero ps and px, and is defined by the properties of the active network 
(i.e. transistor) alone. It should be noted that in general k; takes different values for the different 
parameter environments. If k; < 1 the network is potentially unstable and cannot be conjugate 
matched for MAG; thus in fig. 1, k;j < 1 for frequencies between fa and ff. 


2.2 Transducer gain 


For a given value of k, exceeding unity, the transducer gain G (power output in load termi- 
nation divided by maximum available power from source) may be maximised, so defining 
corresponding optimum terminations ps, py, in terms of k. The terminations are such that the 
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total port parameters (ps + pi), (px + po) are tuned to real values with some degree of mismatch. 
However, the maximised transducer gain G for a given k value is in general different for each 
parameter environment. That this is so may be seen by examining two particular situations of 
practical significance, namely: (a) k—>1 corresponding to a high-gain narrow-band amplifier, 
and (b) k>1 as in amplifiers of very well-defined gain and, usually, wide bandwidth. It may be 
shown that, with little error [3]: 


(a) IF k—1, 


d .\N1/272 
G) Gmax,k & 4 FE, a ee L | 
tie 1)? 


(b) FkS1, 
pai| Î 2E 


(4) Gmax,k © 4 : e 
piz IN ROULE EM) 


A 


Although the quantity | p21/p12 | is an invariant, ki and 2L/(L + M), and therefore Gmax,r 
depend on the parameter environment in both the above cases. 


À more useful parameter than k for the purpose of expressing gain performance is the 
«performance factor» n , defined [3] as 
_(L+M) 


(5) tr) / 
L or 


Thus in the situation of considerable mismatch at the ports, corresponding to case (b) above, the 
transducer gain maximised for à given value of n becomes simply 


sk for k (or n) > 1 
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(6) Gmax,n & 4 Pal 


p12 


and is approximately an invariant with change of the parameter system. In this sense, therefore, 
the performance factor may be regarded as a basic measure of gain capability. 


When unilateralising feedback is applied to the active two-port network in the manner of 
fig. 2, a procedure adopted particularly when high gain is required—requiring the condition 
k—1 in absence of feedback—, the resulting optimised MAG may again be expressed simply in 
terms of the performance factor. Assuming feedback arranged in one of the four environments 
such that the overall network is unilateral (p12 + pr12 — 0), and the MAG is maximised to yield 
the value Goptuniatthen with only small error [3] 


il 
(7) Gopt,unilat UT 
IH; 


assuming | par | > | p12 |. The parameter #; is the «inherent» performance factor of the basic 
active two-port network (without feedback), defined as 


(8) nr 0 11022 
[0 de 
IE 


and is the value of # for zero value terminations. The question as to the best feedback scheme to 
be adopted is readily answered in terms of the above gain expression: Gopt,unilat 1s à maximum 
for the feedback «environment» for which ni, expressed in terms of the corresponding parameter 
set, is a minimum. The table below cites a particular example in which h, z and y environment 
feedback are compared for the common base (C.B.) and common emitter (C.E.) configurations. 
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OC44: 500 kc}s | 


EME Feedback Environment 
ns h 0 0 } 
Configuration hi Gmax,unilat | fi Gmax,unilat | Hi Gmax,unilat 
|.007 404B |.035 3348 | 45 22 dB 
| 41 29 dB 20 32 dB | .029 41 dB 


In this case it is seen that the common-base and common-emitter configurations give similar 
optimum gain performance, but that this corresponds to different feedback systems. 


Equations (6) and (7) may be combined to give 
(9) Gmax,n ‘1 & 16 Goptunilat "1 for n > 


A mismatched active two-port network may thus over the range of performance factor for 
which 


n < 16; 


develop à larger transducer gain than the optimised MAG with unilateralising feedback applied. 
Eqn. (9) summarises compactly the relative merits of simply mismatched, and conjugate matched 
unilateralised amplifiers, in terms of the performance factor. 


3. An invariant stability factor 


The parameters k and n have been seen to be appropriate for defining stability conditions and 
gain performance respectively. In general, however, the values of both factors depend upon the 
chosen parameter environment for given terminations, though the performance factor approaches 
an invariant when n 1 corresponding to large degree of mismatch at the ports. 

k and n may be related to a more fundamental quantity having the property of strict invariance 
[4] in terms of which the results of the previous section may alternatively be expressed. This new 
parameter is accordingly defined as the «invariant stability factor» and is defined in terms of 


(10) n = (2 p11 022 —M)/L 


the value of which may readily be shown to be independent of the chosen set of network 
parameters, and to be between the extremes of —1 and co. The invariance of 1 may be verified 
by writing the MAG of the active two-port network as 


il 
n + M2 —1)12 


It follows therefore that for MAG to be à meaningful expression 


(11) MAG | pe 
| P12 


Neal, 
so defining absolute stability; otherwise if 
—1< ni 


the network is potentially unstable and any transducer gain possible. 
The invariant stability factor is now defined [4] as 


(12) Si = n + (n?—1)1/2 
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and for absolute stability with any, including zero value, terminations 
Si >1 


corresponding to the condition 
k; Mi 


in terms of the inherent Sfern’s stability factor. If the active two-port network is potentially 
unstable, the value of S; is complex, lying on the unit circle on the Argand diagram. 
In general, including the terminations at the two ports, if 


(13) 9 = 2(p1 p2—M)/L 
the network is stable with any values of 6s, 61, if the invariant stability factor 


(14) S= 1 +(n2—1)1/2>1 


: . 
corresponding to Stern's condition 


k>1 
Note that both conditions give 
2 01 02 = L+M 
for Sr 


Thus the invariant factor $S can take the function of Stern’s factor k for the definition of conditions 
of stability. 
As seen above, S; has particular significance in regard to conjugate matched gain performance 


of an absolutely stable system (Si, k; > 1). It may be shown [4] that with conjugate matching for 
MAG 


(15) SE 


where p1, pa are for the values of the terminations concerned at the ports. Thus 


| pai 


4 1 | p21 
| p12 


(16) MAG 
Si | p12 


| Gie| 


where Ge may be defined as the «internal loop gain» under conjugate matched conditions: 
i.e. for any ps, px, the internal loop gain is 


(17) = p21 P12 
p1 P2 


It should be noted, however, that only under the special conditions of conjugate matching is 
| Gil an invariant with change of parameter system. The connection of the factor S; with the 
internal loop gain is not unexpected; the above definition of absolute stability is closely associated 
with the requirement that, in accordance with Nyquisf’s criterion, the locus of Gj in the complex 
plane must not encircle the point (1,0) for any terminations at the ports. 

If the active network is such that either of the inherent factors ki, Sÿ à 1, then 


2 ; 
n > land one Or 2ni, 


611022 > | pie par | 


and for any terminations 
|G| <1 


Under these conditions, Sj — 2n and 
4 , 
(18) Sr RE = 4m 
The effect of the parameter p12 on the total port parameters is then negligible, so that with 


conjugate matching 


e1 S2p11 | 
’ o2 © 2 p22 | 
an 
Î 6 
(19) = mes pee | _ 
e | Gie! | Pi2pal 
giving 
sa 
(20) MAG & par | — where # & 4 nj and ni > 1 
| P12 11 


In the case of a conjugate matched network with negligibly small reverse signal transfer, there- 
fore, the invariant (inherent) stability factor S; may be identified with the performance factor 
n & 4 hi. 

In general, with any values of ni, Si, but with terminations at the ports such that 


kn 1 


Ho 
(21) SR = 4n 


In the situation of large k or n, the performance factor may thus be interpreted simply in terms 


of the invariant factor S. If, further, conditions of considerable mismatch exist at the ports, 
corresponding to # à 1j, the maximised transducer gain for given large S is, from eqn. (6), 


Î 
+ — for S SA 


| 
| S 


(22) Gmax,S Æ 16 


| p12 


4. Instability and oscillation 


In each of the four basic feedback systems of fig. 2 the passive feedback network may be 
adjusted for «maximum instability» in the sense that the largest value of the real part o1, of the 
load termination for which the overall network is unstable is maximised. The active two-port 
network may then be unstable and oscillate with the greatest possible output port loading. The 
conditions for maximum instability may be determined by examination of the total output port 
parameter i.e. output mesh impedance, or admittance across the output terminal pair, as appro- 
priate: the real part of this parameter, negative for instability, is required to become zero for the 
largest possible value of pr. 

Consider the y-feedback environment shown in fig. 4, in which the source admittance 
ps = ys — 0. The total output port parameter pyr is the sum of the ouput admittance yo and load 
admittance y. For instability and oscillation, the output conductance go needs to be negative 
with 


(23) (— go) = gx > 81 
while the frequency of oscillation is defined by the condition 
(24) bo = — br, 


for zero output port susceptance. 
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Fig. 4. Application of feedback to cause instability 


(a) y-environment (shunt-shunt feedback): for instability 
go + g5 < 0: 
(b) Two forms of feedback network 


(a) 

ñ Til 2 SE 2 

© [jB+ | = 

TRE Ù Un 
@) 


Thus, maximum instability is attained when the negative output conductance gx is a maximum. 

Fig. 4 shows two alternative ideally lossless feedback networks, the circuit elements of which 
may be determined for maximum gx. In terms of the transformer feedback network, for 
maximum instability 


DEN 
dT 
(25) 
LEE 0 
dBr 


So defining optimised feedback network parameters T°, Br° and corresponding maximum 
negative output conductance gx °. The values of these optimised quantities are [5]: 


g21 — £12 
Br° = —b11 en | —- . 


bo == b12 


| ne jun 
(26) D Del) | 
2 [g11 (g21 — g12) + b11 (b21 — b2)] 


il . 
and gN° = ——[(go1 + g19)2 + (bar — b12)?] — go 
gi1 


In the above expressions, the lower-case symbols g11, b11, etc all refer to the y-parameters of the 
basic active two-port network. Similar results may be obtained in terms of the appropriate 
parameter sets by analysis of the other three feedback systems of fig. 2. 

Alternatively, the feedback network may be a pi-configuration of susceptances, as shown in 
fig. 4, and this is often more convenient than the transformer network in practice. The suscep- 
tance components are readily evaluated from the relationships [5,6]: 


Bi Dr) 
(27) Ba — br TPBr- D 
B3 = Br T 


Many other forms of feedback network are, of course, possible, and may be derived from the 
results given above. 
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In the case of a junction transistor, the form of the feedback network for maximum instability 
is sensitively dependent upon the phase angle of its forward transfer parameter (y21 in the above 
analysis). Thus if the transistor is in the common-base configuration as in fig. 5 the phase lag of 
the short-circuit current gain parameter « of the transistor or ofits internal current gain parameter 
a” in the equivalent circuit shown in fig. 6, is most significant. Writing a” as 


Aie JV di) Co 

it is interesting to see how the nature of the optimum feedback network changes as the phase lag 
D increases; it should be noted that under practical conditions of oscillation DA, can be large, 
possibly lying between x and 3x/2 at the maximum frequency of oscillation. The forms of the 


ES DITION TRANSFORMER FORM 


o<d, <4 
; 
Res 1) 
4-4 
1 (=) 
= Pr LA 7 
T>1) 
Internal current qœin a’ 
of transistor, defined in 
Fac. FA 2 
a'= pp -;Ÿ (T =) 
CA > 7 
TJ AE 
Ë [ (T: -ve) 
(4 B, : incuded in Ca. 
-Ÿ 
then jatl? = 2u 


de = + 


1 


Fig. 5. Dependence of feedback network forms on the phase of transistor gain: common base configuration 


transformer and pi feedback networks are tabulated in fig. 5 for various phase lag conditions 
(see [6]), showing how in general these networks change as the frequency of oscillation increases. 
Two conditions are of particular interest: for 


Py =D, where|a|2=2u 
T = 1 and the transformer, and one element of the pi-network, are redundant; for 
Da == r/2 


T = © and the transformer feedback network cannot be realised, though the element 
values of the pi-network are all finite. 


4.1. Maximum frequency of oscillation 


The maximum oscillation frequency fmax may by derived from the expression of eqn. (26) 
for maximum negative output conductance gx by invoking the condition 


(28) [ gn = 0 
at | = fmax 
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In this respect it should be noted that the feedback network has in the foregoing discussion been 
assumed lossless, as required for oscillation at frequencies up to finax. Thus, solving the equation 


(29) (go1 + g12)2 + (bai — b19)2 — 4 g11 ga = 0 


in terms of frequency, where the y parameters refer to the transistor in any configuration, 
yields fmax. Expressions for the y parameters of the transistor may be derived from the equivalent 
circuit representation of fig. 6, in which the most significant elements are the internal current 


Fig. 6. Equivalent circuit of the transistor at high frequencies. 
generator aie and the extrinsic components rpp’, Ce. If the assumption may be made that the 
emitter diode resistance ra is negligibly small, eqn. (29) reduces to the equality [6] 
| a" |? 
TV lb! Ce 


at  f— fmax 


fm may be determined by finding a point on the locus of 4’ at which the above condition is 
satisfied at the frequency concerned. The assumption of suficiently small rq in the derivation of 
eqn. (30) implies that simultaneously [6] 


(31) ra < (1— 4) mp [1 + (ora Ce)?] 


{4 
V Thb 


(30) Îm = 


À 


ul - [1 + (w ra Ce)?] 


oCe 
vrpp' [1 + (wo ra Ce)?[? 
v + (1 + pu) w ra Ce 


where Ce = Cte + Ce. It is evident that these assumptions and eqn. (30) break down for v — 0. 


If it may further be assumed that 
1 


Oo 

1+jflfa 
a form of frequency dependance that is however rarely obeyed by a transistor, particularly a 
device of drift type, then explicitely, 


/ 


( ) fm 8 T lbb CS 
Although this last expression is frequently quoted as an indication of the maximum frequency of 


oscillation, it is clearly inaccurate in general, possibly grossly so, owing to the excess phase lag 


exhibited by a’ at f = fmax. 


4.2. Power output 


The optimisation process discussed above ensures that oscillation persists for values of output 
load px, up to the maximum possible, at a given frequency. If the output susceptance of the 
transistor with feedback applied is negligible compared with that of either the capacitance or 


13 


inductance elements of a frequency defining tuned circuit at the output port, optimisation of 
feedback for maximum negative output conductance also results in the maximum rate of build 
up of oscillations — as may be required, for example, in a gated oscillator or super-regenerator. 
Such a feedback network does not, however, allow efficient operation of an oscillator from the 
point of view of power output developed in a load of less than the maximum possible value, for 
oscillations are found to limit at very small amplitude. 

It is possible alternatively to optimise the feedback network so that maximum output power 
is developed with approximately class-A operation. Limiting of oscillations must be arranged to 
occur simultaneously at the emitter and collector diodes, à situation that may be considered to 
exist when the peak excursions of emitter current and collector voltage are equal to the mean 
emitter current Ze and collector supply voltage V4 respectively. The collector efficiency then 
approaches 50%. Assuming that such limiting action occurs and that the circuit waveforms are 
approximately sinusoidal, corresponding to the transistor exhibiting the small signal parameters 
of the mean operating point (Je, Ve) over the whole oscillation cycle, the load power may be 
evaluated and maximised by best choice of the feedback network. For the transformer feedback 
arrangement with the transistor in the common base configuration of fig. 4, the optimised net- 
work parameters are [7]: 


(33) “ 


= É (gu/ly11|?) | T° 
— Ve la 
BE ni 
where T°, Br° are the parameter values for maximum negative conductance given by equations 
(26), y11 refers to the transistor alone as before, and & is the complex emitter-collector current 
gain. The quantity g11/|y11l2 — R (1/y11) — R (h11) is the real part of the input port (emitter- 


base) short circuit impedance. The corresponding power output and load conductance values 
are [7]: 


noie Ie? 811 
PL - | 
2 2001|ÿ11|2 
4 | Ze (2 AI 2 
(34) ne 2 al —| és si = | e (g11/| yuil . Pa 
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Thus in order to realise the maximum class-A power output, only the transformer ratio of the 
feedback network requires to be changed from the previous optimised value. The load conduc- 
tance, given by eqn. (34), is much less than gx ° in eqn. (26). 


4.3. Frequency stability 

Design of à transistor oscillator for good frequency stability involves considerations quite 
different to those discussed above. The frequency defining condition of eqn. (24) needs to be made 
to a high degree independent of changes of transistor parameters (viz. changes due to supply 
potential or temperature fluctuations, device ageing etc). Thus the rate of change of load sus- 
ceptance with frequency, dbr/do should be made as large as possible, requiring a tuned circuit 
of large C/L ratio; also possible change in the output susceptance bo of the transistor, with feed- 
back applied, should be minimised. Now the output susceptance for the feedback arrangement 
of fig.4, with the transistor in common-base configuration, may be written [7]: 


KE LME vT 


35 bo = I(h22 Eine ee: LISE 
À : PUR RDA + KO RD + 9 
Where 

(36) ke Lin) —1/Br 


R(h11) 
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It is convenient here to use h-parameters for the transistor, so that 99 is the open-circuit collector 
admittance. The first term on the right hand side of eqn. (35) thus depends only on the collector 
capacitance C& for the common base configuration; this term would be in general much larger 
and more subject to change for the common-emitter configuration. The second term may be 
made zero if Br is chosen to make K zero at the frequency of oscillation. This term becomes non- 
zero, however, if I (h11) subsequently changes and then depends also on R(h11). The effect of 
R(h11) may be made adequately small by adding a resistor of sufficiently large value in series 
with Bp; the magnitude of the second term in eqn. (35) is then reduced, as is also the maximum 
frequency of oscillation of the transistor. The value of I(h11) falis as the ratio of the frequency of 
oscillation to the current gain cut-off frequency f, is reduced. The third term of eqn. (35) depends 
on y, 1.e. the imaginary part of «, and this also is therefore reduced with decrease in the oscil- 
lation frequency. 

Application of the above measures, in particular adoption of large load C/L, low frequency of 
oscillation and addition of resistance in the feedback path, lead to good definition of frequency, 
having initially chosen common-base as the most suitable configuration of the transistor. 


5. In conclusion 


Considerations relevant to the application of the junction transistor as an active two-port 
network have been discussed, with emphasis on fundamental rather than detailed issues. Investi- 
gation of the conditions for stability or instability of the device with arbitrary terminations at the 
ports has led to the expression of transducer gain performance of amplifier configurations, and 
optimisation of oscillators in regard to regeneration, power output or frequency stability. 

It is worth noting that most of the results obtained are quite basic, applicable to any active 
two-port network ; the transistor may be regarded as a particular example of such a network, the 
study of which has resulted in conclusions of wide significance. 
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Optimization Procedures for Sampled-Data and 


Digital Control Systems 


by E. I. Jury 


Paper presented at the 5th SFAC Meeting (Swiss Federation of Automatic Control) 
in Lausanne on May 20, 1959. 


Optimisierungsmethode für «sampled-data»- und digitale Regelsysteme. Im Hinblick auf ihre grosse Anpas- 
sungsfähigkeit in Bezug auf die jeweils benGtigten Übertragungseigenschaften gewinnt die Verwendung 
von Digitalrechnern in Regelsystemen ständig an Wichtigkeit. Die vorliegende Arbeit behandelt die 
Hauptzüge der Optimisierungsmethode für die Anwendung auf solche Systeme. Ein typisches Block- 
schaltbild für ein solches System ist in Fig. 1 gegeben. Das Optimisierungsverfahren beruht auf der 
Methode der kleinsten Fehlerquadrate, wobei der Fehler als Differenz zwischen vorhandener und ge- 
wünschter Ausgangsgrüsse definiert ist. Während für rein diskrete Systeme die z-Transformation ver- 
wendet wird, kommit für gemischte (diskret-stetige, z. B. ein Digitalrechner in einem stetigen System) 
die modifizierte z-Transformation als Verallgemeinerung der z-Transformation zur Anwendung. Dabeï 
werden natürlich für die Berechnung der digitalen Übertragungsfunktion D*(2) die Bedingungen für 
physikalische Realisierbarkeit und Stabilität im Optimisierungsverfahren berücksichtigt. Einige der wich- 
tigsten Eigenschaften der modifizierten z-Transformation sind aufgeführt. Zur Veranschaulichung der 
erhaltenen mathematischen Ausdrücke wird in einem Beispiel die Anwendung der Methode auf einen 
praktischen Fall gezeigt. Schliesslich gibt der Schluss der Arbeit Angaben über die Ausdehnung der be- 
trachteten Optimisierungsmethode für die Behandlung von Regelsystemen mit Digitalrechnern als 
Kompensations- und Filterglieder in Rückkopplungszweigen sowie auf andere Systeme (z. B. solche mit 
endlicher Impulsbreite). 


Procédés d'optimisation pour des systèmes «sampled-data» et digitaux de réglage. Grâce à ses grandes possibilités 
pour le choix des propriétés de transfert, l’emploi du calculateur digital dans des systèmes de réglage a 
gagné beaucoup en importance. Le travail présent s'occupe des principes de la méthode d’optimisation 
pour l’application à ces systèmes, pour lesquels la fig. 1 donne un schéma de principe courant. Le procédé 
d'optimisation se base sur la méthode du carré moyen minimum d’erreur, cette erreur étant définie comme 
la différence entre les signaux de sortie actuel et désiré. La transformation en z (transformation de Laplace 
adaptée aux systèmes «sampled-data» par l’introduction d’une nouvelle variable z — eTS), appliquée aux 
systèmes purement discrets, doit être généralisée pour l’application aux systèmes mixtes en question (sys- 
tèmes discret-continus, par ex. un calculateur digital incorporé dans un système continu): on obtient la 
transformation de z modifiée, dont quelques propriétés importantes sont mentionnées. Dans le calcul de 
Ja fonction de transfert digitale, 1l doit évidemment être tenu compte des conditions de réalisation physique 
et de stabilité. Les relations mathématiques dérivées pour le comportement des systèmes mixtes sont 
illustrées à l’aide d’un exemple typique. Le travail se termine par des indications sur l’extension de la 
méthode en vue de son application à des systèmes comprenant un calculateur digital comme élément de 
compensation ou de filtrage dans une boucle de contre-réaction, ainsi qu’à d’autres systèmes. 


The use of digital computers in control systems is becoming increasingly important in view ofits flexibility 
in types of responses and filtering obtained. This paper discusses the main features of optimization proce- 
dures for such systems. À typical block diagram of the actual system, the desired system and the system 
error is indicated in Fig. 1. The optimization procedure is based on minimizing the mean squared error 
(difference between actual and desired output). Physical realizability and stability conditions are included 
in the optimization procedure for obtaining the digital transfer function D* (2). The method of approach 
is based on the modified z-transform and certain features of this method are briefly indicated. An example 
is included to illustrate the various mathematical terms discussed and to indicate the application of the 
method to a typical system. The extension of this method to the use of digital computers as means of 
compensation and filtering in a feedback path and to systems with pure delays and of finite pulse width 
are indicated in the conclusion. 
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1. Introduction. 


In recent years the Wiener Theory of prediction and optimization has been applied to sampled- 
data or digital control systems. Similar to the continuous case the approach is based on the 
optimization of the error (difference between actual output and desired output), when the 
correlation function or the power spectral density of the input are given. For strictly discrete 
systems (systems in which response at sampling instants is only required) the z-transform 
theory is applied. 

The z-transform theory is widely developed for the analysis and design of linear sampled- 
data systems. If the Laplace transform of à continuous time function f (f) is denoted by F (s) = 
L{f(#)], then the z-transform of F(s) denoted by F*(2) is given by the following integral 
equation : 

i c+joo 1 
F# (2) = Z [F9] = Poe d 


2x; RE) 


CCS 
In some cases F(s) has only simple poles, then its z-transform can be written as: 
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(s) BO 
SSD Sa ON are the simple roots of B(s) — 0 
Pi LE ;  T'= sampling period. 

ds |s=5n 


For other forms of F(s), tables of z-transform [3] can be used or the above integral can be 
similarly evaluated. 

The theory to be applied for the case of mixed discrete-continuous systems (when a digital 
computer is incorporated in a continuous system) is the modified z-transform method. The 
modified z-transform is an extension of the z-transform method for obtaining the system 
behaviour at all instants of time. 

Similar to the z-transform, this method can be written in the form of an integral which 
yields the modified z-transform of F(s) and denoted as F*(7, M) — Zim [F(s)|. It is given by the 
following relation: 


c+joo . 
Mn 7 [Pole i ro ae a 
MT 2x j Le TR PR 


When F(s) — A(s)/B(s) has only simple poles, the evaluation of the above integral vields: 


A(sn) eMsnT | 


Lt Ne 


FX (z,m) = z 1 _. 0<m<l 

General tables of modified z-transform [3] are available which can be easily used for more 
general cases of F(s). 

In this paper the modified z-transform approach is applied and extended to mixed systems, 
whereby the optimization criterion based on the minimization of the mean square error 
(difference between actual output and desired output at all instants of time) is applied. 
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The technique is discussed only for the feedback configuration shown in Fig. 1, however its 
extension to other configuration is straightforward. 


n(t) RÉ 
one AR Te tu. 
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Digital 
Compensator 
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Desired 
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Theoretical 
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Fig. 1 Block diagram of the actual system, the desired system and the system error. 


An example is included to illustrate the application of the theory to a particular system and 
to indicate the various statistical relationship required and the important assumptions under- 
lying the optimization procedure. 


2. Optimization Procedure. 


We will confine ourselves in this discussion to the optimization procedure for series discrete 
compensator in the forward path as shown in Fig. 1. 

Extension to other configurations can be easily applied. For the system in Fig. 1, the given 
r(t) and n(f) 1. e. the input signal and the noise are assumed stationary and random. The desired 
transfer function is L(s). The overall transfer function K(s), which contains the plant G(s) and 
DX*(2) discrete compensator is determined to obtain the minimum mean squared error and to 
yield a stable feedback system. The latter condition requires that K*(z,m) = Zm[K(s)] should 
contain all the zeroes [1] which are outside the unit circle of G*(z,m) = Zm[G(s)]. 

The actual output c(#) is the response of the system for the input Ri(s) = Rs) + N(s), and 
the desired output ca(t) is the response for the input Rs) to the system defined by the transfer 
function L(s). The difference between the desired ouput and the actual output is called the 
error, and the system is to be designed for the minimization of the mean squared error quantity. 


Mathematically the error transform is given by: 


(1) HO 10 CN 
where 

@) Cat) LE ROUE 7: 
and 


The actual output is given by: 
(3) C(s) = K(s) R1* (2) = K(s)[R*(2) + N*()], 2z=e7s 
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Taking the modified z-transform of equations (2) and (3), we obtain: 


(4) Ca* (2, m) = L*(2, m) R*(2) 
(5) CX(2, m) = K*(2, m) [R*(e) + N*(2)] 

From the physical realizability condition, K*(z, m) can be written as: 
(6) K(2,m) = ZE P* (2, M) F*(2) 


where 2-1 P*(z, m) is the numerator of G*(z, m) = Zm[G(s)] and F*(2) is an arbitrary func- 
tion which is to be determined from the minimization procedure and has no poles outside the 
unit circle. 


3. Mathematical theory and optimizations. 


Having defined the optimization procedure in the preceding, we now apply the mathe- 
matical theory of the modified z-transform to obtain F*(2) satisfying the above conditions. 
This can be done by starting with the modified z-transform of Eq. (1) and using Eq. (4) and (5): 


(7) EX(z, m) = Ca*(, m) — CY(e, m) 
= m ZIP mn) FORGE PPXC MIEAOINRE) 
The power spectral density of the error e(f) is defined by: 
(8) DEC m) EEE MIE EL #)] = 
k = 00 1 q=N 
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The notation E is used here as a symbol operator to indicate the limiting process of the last 
term in the above equation.The mean square error is given [1] by the constant term of D6*(z, m), 
averaged over m, that is: 


—< 1 ————— dz 
(9) A D. Ÿ DEN") PR 


unit circle 
Therefore our principle objective is to find a function of z, F*(2) that will minimize the 
(10) constant term of (Dee* (z m)). 


Since by definition E [R*(2) R*(z-1)]1) = D,r*(2) (the power spectral density of the 
sampled input), therefore equation (8), when equation (7) is substituted for E*(2) and using 
this definition, yields: 


Dee*(z, m) = [L*(&, m) — 2 1PX (2, m) FX (2)] [L* (GEL, mm) — 2P* (LL, m) FX (2 D)] Drr* (2) 
(11) PC mP TEE) PCA), 
— 2 [L* (2, m)—2Z1P* (2, m) EX (2)] PX (2 LL, m) FX (271) Dar * (2) 


— 21 [L* (2, m) —2P* (21, m) EX (2 2)] P* (2, m) FX (2) Drn* (2). 
1) E[R*(2) R*(2-1)] = ®;r*(2) = lim SRI Rx*(2) Rw*(z71) 


N—00 
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Assume F*(z) as the sum of two functions: 
(12) FX (2) = Fo*(z) + F1*(2) 


where Fo* (2) is the optimum function which gives the lowest constant term of (Dee*(z, m)) 
and F*(2) is any function included to satisfy the realizability condition. F1* (2) is the difference 
and also physically realizable (no poles outside the unit circle). Eq. (11) can be written in terms 
Of Fo* (2) and F1*(2) as: 


(13) Dee*(z, M) = Deeo* (2, M) + Po*(e, m) + Wi1*(e, m) + Fo*(z, m) 
where 
Deeo* (7, m) = the right-hand side of Eq. (11) with F*(2) 
and F*(21) replaced by Fo*(2) and Fo* (21) respectively 
(14) Vo*(z,m) = PX (2, m) P* (21, m) F1* (2) F1* (22) [Drr* (2) + Pnn* (2) + Prn* (2) + 
+ Dnr*(2)] 
(5) Wi*(e,m) = [21 PX (2, m) Fo* (2) —L* (2, M)] [Drr* (2) + Dar* (2) P* (TE, m1) F1 * 
1) + P* (2m) P* (TL, m) Fo* (à) F1* (21) [Drn”* (&) + Pnn*(2)] 
(16)  Fa*(e, m) = [P* (21, m) Fo* (21) —LX (EE, m)] [Drr* (2) + Prn* (2)] 271 P* (2, m) 
F2 (2) LP M) PE mn) CE) 70) [POP AAC) 
It is noted from Eq. (14) that Wo*(z, m) is the auto power spectral density of P*(z, m) 
F1*(2) [R*(2)], and therefore the c. t. {Wo*(z, m)) is positive, unless F1*(2) = 0 (a trivial case). 
Therefore, the necessary and suficient condition for 


(17) ct. (Dee* (z,m}> ct. (Deco*(z, m)) 


with arbitrary F1*(2) is that 


1 
(18) 1! [et (PiX (2, mm) + 'o* (2, m1) )] dm — 0 
0 
Comparing Wi*(z, m) and Wo*(z, m), the following symmetrical property can be written: 
(19) BAC CL in) 
Thus 
(20) ce mi=celons n)] 


and therefore, Eq. (18) reduces to 
1 
(21) rl [c. t. 1X (2, m)}] dm = 0 
0 


It can be shown [1] that the constant term of a certain function of «z» which can be expanded 
in power series is equivalent to the residue at the origin of that function divided by «2» provided 
the contour of integration is taken as the unit circle around the origin. Thus: 
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or from Eq. (21): 
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1 
1 dz IR 
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Substituting for W*(z, m) from Eq. (15), the integrand of Eq. (23) becomes: 


(24) Z LP * (2, m) = F1* (21) (PAG 1m) PL ml Ce) nn re) 
— [LX (2, m) PX (21, m)] [Drr* (2) + Pnr* (2)]) 


where 
(25) Drir1”* (2) … Drr* (2) + Dan* (2) 7 Dire (@): 


It should be noted that in Eq. (24), F1* (2-1) has no poles inside the unit circle, since F1*(2) 
(from physical realizability) is assumed to have no poles outside the unit circle. 


We introduce a new function X*(z) as follows from Eq. (24). 
(26) HART E) SENTE) 


If X*(2) has no poles inside the unit circle, then Eq. (23) is satisfied. This is also a necessary 
condition due to arbitrariness of F1*(2). If X* (2) has any poles inside the unit circle, a function 
F1*(2) ca be selected such as the sum of the residues in Eq. (23) will not vanish and the optimum 
condition will not be satisfied. 


From Eq. (24) and (26), X*(2) is 


(27) KA) PAC MAPS PER M) Drirs () Fo) 74 
— LX(z, m) P*X(2 1, m) [Or* (2) + Dar* (2)] 


Due to symmetrical property of D;,r,* (2), it can be factored into the product of Y*(2) 
and its conjugate as: 


(28) Drirt(e) = V* (2) Y* (2-1) 
where 
q : 
+ (1— Bi 71) 
(29) VOA 


[le] 


(1 — Gi 2) 


(30) CRE 


LAS 


Au: 2) 


em 
en 


Therefore Y* (2) Y*(2-1) represents the product of factors of Dir. * (2) with zeros Gi or 
poles 3 inside and outside the unit circle respectively. 
We can perform the following factorization from the product of the modified z-trans- 


form [1]: 
(31) PX(z,m) P*(2"1,m) = B*(2) B*(z-1) 
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where B* (2) and B* (2-1) have all its zeros and poles outside or inside the unit circle respectively. 
If we let 


(32) TE MmPEE 1 10 
then Eq. (27) can be simplified using relations (28), (31) and (32): 
(G3) X*()=B* C1) V6) EDEN C) AB CM) PER EEE SC) 
Dividing Eq. (33) by B*(2) Y* (71) we obtain: 
TX (2) [Drr* (2) + Dnr* (2)] X*(2) 


1B*(>-1) Y* À L i 
(34) z-1B*(>1) V*(2) Fo*(2) B*(2) Y* (71) BOY 


Let Q*(2) represent the following: 


TX (2) [Drr* (2) + Pnr*(2)] 
35 X (2) = — . — 
(35) Q*(e) Bé) v4( 2) 


In general, Q*(2) can be written as a partial fraction expansion in terms of roots À; of the 
denominator and à polynominal II* (2). Further it could be divided into two parts depending 
on whether À; is inside the unit circle or outside of it, as follows: 


(36) OX) = O1) 00" E) 
where 
(37) Qi* (2) = _ 

Per 
38 * (2) — al 
(38) Qo* (2) 2 nn (2) 


Substituting the above for Q*(2) in Eq. (35) and inserting in Eq. (34) yields: 


(39) a DO Di Ou ie 
B* (2) Y*(21) 


It is noted that all the poles in the left-hand side are inside the unit circle while the poles of 
the right-hand side are outside the unit circle. In order to satisfy Eq. (39) both sides of this 
equation should be zero, taking into consideration the fact that Q;*(2) has no constant term 
and Fo*(2) has no z-term. Thus Fo* (2) is obtained from the following equation: 


(0) 1 BA (e71) VA (2) Fo* (5) — Q*(à) = 0 
OT 

x 
qu) Fox (e) = 20 @) 


B* (271) Y*(e) 


This is the optimum F*(2) [denoted as F5*(z)] which minimizes the mean squared 
error e(t)2, the optimum design criterion used in this paper. 


To obtain the optimum digital filter (or compensator) Do* (2), we find the optimum K*(z, m) 
denoted Ko*(z, m) as follows: From Eq. (6) and (41): 


Primo 


(42) SE. CASA 
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From Fig. 1, the total transfer function is given: 
DX(2) G*(z, m) 
1 + D*(2)  G*(2) 


(43) l'ENA 


Thus, the optimum compensator Do*(2) in the forward path is given by: 
Ko* (+, m1) 


(44) Do — | 
di GX*(2, m) — GX (2) Ko*(z, m) 


where Ko*(z, m) is given in Eq. (42). 
There exist several methods for the realization of the discrete compensator [3]. For brevity 
of the discussion, these are not included in this paper. 


4. Extension of the optimization theory. 


The preceding optimization techniques, although applied to a specific feedback system, can 
be extended to other configurations in particular to the system shown in Fig. 2. 


n(t) K* x 
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Fig. 2 Control system with digital compensator in feedback path. 


Further extension of the techniques can be applied to plant transfer function with pure delay, 
to multiple sampled systems and finally to sampled-data systems with finite pulse-width. 
Although, each of these cases is to be treated separately, however the optimization criterion is 
the same. 

Of importance is the comparison between digital filtering (using discrete elements) or 
analogue filtering (using linear networks) and the advantages of each method should be in- 
vestigated. 

Extension of the technique to non-linear or to time-varying systems should not be over- 
looked. Indeed, this extension offers many advantages although the theory becomes much 
more involved. 


5. Illustrative Example. 
To illustrate the application of the optimization techniques developed so far, the method 


is applied to the system shown in Fig. 1 with the following constants: 


1 
Fa Pre) = PET 
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(46) Dan*(s) = 1 
Drn(s) = Pnr(s) — 0 
The transfer function of the system (unalterable component)?) 
47) CO = TE 
ALT 
Assume the desired transfer function of the whole system as a pure integrator (it can be a 


predictor e*8 or unity or any other desirable transfer function): 


(48) L(s)=—, K=— constant 


From Eq. (6), the overall transfer function should be of the following form: 
(49) KT HP ES 
where P*(z, m) is the numerator of z + G*(z, m) and F*(2) is an arbitrary function of z which 
should have no poles outside the unit circle (physical realization). 

GX*(z, m) is obtained from G(s) using tables or the convolution integral?) to yield: 
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(50) Cm) — 


(Steel 
therefore the numerator is 
(51) PACE mn) er 


From Eq. (31): 
1 


(52) Be) — 1e m) dm = 0,632 
0 
The modified z-transform of L(s) is 
Kz1 
(53) TE 
1— 21 
Then T*(2) from Eq. (32) equals 
0,632Kz 1 
(54) T*(2) = ; en 


To find Y*(2), obtain z-transform of Eq. (28) 
1 


EE) 


2?) Usually in a practical system a zero-order hold circuit is incorporated in G(s) 


(55) Drr* (2) = 
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Since ®,r*(2) = 1, thus, from Egq. (25) 


1 1—(W2—1)21 1—(/2—92 
D— V2 1-1) 112 
Factoring the zeros and poles inside the unit circle: 

1 — 0,414 3-1 


(57) VX (2) = 1,300 
10707 


(56) Prim* (2) — 
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QX(2) is obtained from Eq. (35): 


bn cree) à | O368k71 
B*(2) Y*(7 1) (1— 3-1) (1 — 0,707 2-1) (1 — 0,414 1) 


Qi*(z) from Eq. (37) consists of these terms of the partial fraction expansion of Q*(2) 
which has poles inside the unit circle. Therefore, 
2 = tt 2 le. 
(59) Qr* (2) = 2,120 Kz 1,256 K- _ 0,864 (= 02820 KR 
1— 231 1 — 0,707 2-1 (1— 2-1 (1— 0,707 2-1) 


Thus Fo*(2) from Eq. (41) is 


2Qi* (à) _ K(1—0,282z1) 


60 Ho = — _ 
‘4 2 B*G1)V*() (1—271)4—041471) 


The optimum transfer function Ko*(z, m) is obtained from Eq. (42) 


NES eee 
21) Ko*(z, m) = 2 1P*(2, m) Fo* (2) = Ke (1—0,282z1)e# 
(1=2"1) (10,142) 


The optimum digital filter (or compensator) is obtained from Eq. (44). 
K (1 — 0,282 z-1) (1 — 0,368 7-1) 


(62) Do* (2) 
1—K— (1,414 —0,282K) 1 + 0,4147 2 


For example if K = 0,1, then 
0,110 (1 — 0,282 z-1) (1 — 0,368 2-1) 


(63) Do*(&) = 
1—1,542 1 + 0,460 22 


The step response for such an optimum system is shown in Fig. 3. 


6. Conclusions 


The statistical theory that has been developed extensively in the last decade is applied to 
discrete or mixed systems. The optimization criterion is based on the minimization of the mean- 
square error. It was assumed that the system is linear and time-invariant. By the use of the modi- 
fied z-transform theory the optimization techniques have been developed. Necdless to say that 
other techniques based on the Laplace transform theory and frequency methods could also be 
applied with some advantages. This has been done by other authors [2, 4]. The time-domain 
procedure indicated in this paper is just an extension of other methods. It summarizes the modi- 
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fied z-transform and extends its applications to the actual optimization procedure. Extension 
of these techniques to other pulsed systems is amendable and certainly warrant further investi- 
gations. 
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Die Kreiszylinderdiode im Raumladungsgebiet bei beliebigen 
Emissionsstrôomen und Austrittsgeschwindigkeiten 


Von Peter A. Tschopp 


Mitteilung aus dem Institut für Hôhere Elektrotechnik der Eidgenüssischen Technischen Hochschule 
in Zürich. 


Die vorliegende Arbeit beschreibt zwei numerische Verfahren zur Lôsung der Raumladungsgleichung 
der Zylinderdiode bei beliebigen Enmissionsstrômen und -geschwindigkeiten. Nach einleitender Über- 
sicht über die bisherigen Arbeiten folgt eine kurze Behandlung der Relaxationsrechnung für dieses 
Problem. Der zweite Teil behandelt die elektrische Analogie in Form einer Widerstandskette mit zu- 
sätzlichen Stromquellen. Für ein gegebenes Verhältnis von Anoden- zu Kathodenradius wird damit 
eine vollständige Kurvenschar inklusive Grenzkurven bestimmt für Ort und Betrag des Spannungsmini- 
mums in Abhängigkeit von der normierten Emissionsstromdichte x. 


La diode cylindrique à charge d'espace pour des courants et des vitesses d'émission arbitraires. L’article décrit deux 
procédés numériques pour la solution de l’équation de Poisson relative à ce problème. Après un aperçu 
des travaux antérieurs sur ces questions, l’application du calcul de relaxation à la diode cylindrique est 
brièvement traitée. La deuxième partie de l’article décrit l’analogie électrique formée par une chaîne de 
résistances avec des sources de courant supplémentaires. Pour un rapport fixe entre les rayons d’anode et 
de cathode, un ensemble complet de caractéristiques est donné pour la valeur et le lieu du minimum de 
potentiel en fonction de la valeur « normalisée du courant d'émission. 


The cylindrical space-charge-diode for any initial velocity and current. This paper deals with two numerical 
solutions for the space-charge equation of the cylindrical diode. After a survey on previous publications 
on this subject, the application of the relaxation method is treated briefly. The second part of the paper 
describes an electric analogue computor consisting of a resistance chain with additional current sources. 
For a given ratio of anode to cathode radius, a complete set of characteristics including limiting curves 
is given for the value and for the position of the potential minimum as dependent on the normalized 
emission current density. 


1. Einleitung 


Während die Theorie des stationären Verhaltens der parallelebenen Diode im Raumladungs- 
gebiet nach den grundlegenden Arbeiten von W. Kleen, G. Plato, H. Rothe, M. J. O. Strutt, 
À. van der Ziel u. a. schon früh als abgeschlossen gelten konnte, ist dies für die Kreiszylinder- 
diode heute noch nicht der Fall. Ihren Ausgang nahm die Theorie der kreiszylindrischen Diode 
mit einer Verôffentlichung von 1. Langmuir (Schrifttum [1, 2]), in welcher er zum erstenmal 
die seither als Langmuirsche Beziehung bekannte Gleichung für das Verhalten der zylindrischen 
Raumladungsdiode bei verschwindender Emissionsgeschwindigkeit der Elektronen ableitete. 
Die Langmuirsche Korrekturfunktion 82 war damals noch recht ungenau gegeben. Spätere Ver- 
ôffentlichungen von W. Schottky [3, 4, 5, 6] brachten Erweiterungen und Ergänzungen für die 
Kreiszylinderdiode auf theoretischer und experimenteller Seite, so auch die Berechnung des 
Emissionsstromes bei Maxwellscher Verteilung der Emissionsgeschwindigkeiten für die Elek- 
tronen. Die damals noch «neue» Existenz einer eigentlichen Elektronenemission wurde zur 
gleichen Zeit durch die experimentellen und theoretischen Arbeiten von S. Dushman [7] u. a. 
näher begründet, während spätere Verôffentlichungen wiederum von I. Langmuir und Mit- 
arbeitern [8, 9] die genaue Berechnung der f2-Funktion für den Fall der Innenkathode und der 


Aussenkathode brachten. 
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Erst 1939 beschäftigten sich J. Crank und Mitarbeiter [12] mit dem Fall der Kreïszylinder- 
diode mit endlicher konstanter Austrittsgeschwindigkeit der Elektronen. Es gelang ihnen, mit- 
tels mechanischer Integration mit einem Differentialanalysator für bestimmte Bedingungen 
Lôsungskurven für diesen Fall zu erhalten. Ihre Anwendung erwies sich jedoch als schwierig, da 
zwischen den physikalischen Randwerten und den entsprechenden Kurvenparametern keine 
cinfachen Bezichungen herrschen. 

Nach dem 2. Weltkriege wurde von verschiedenen Seiten an diesen Problemen weiter- 
gearbeitet. L. Page und N. J. Adams [14] gaben 1945 eine weitere Lôsung analog jener I. Lang- 
muirs [1,9] mit einer Potenzreihe für B2, die für gegen Eins grosse Verhältnisse von Anoden- 
radius zu Kathodenradius bedeutend schneller konvergierte. A. van der Ziel Iôste darauf den Fall 
konstanter Austrittsgeschwindigkeit bei Vorliegen eines Potentialminimums mit dem Potential 
Null und schliesslich verôffentlichten L. Page und N. J. Adams [17] eine vollstindige Lôsung der 
Raumladungsgleichung für beliebige Emissionsstrôme und -geschwindigkeiten. Abgesehen von 
der etwas mühsamen und in weiten Gebieten sehr ungenau ausführbaren graphischen Aus- 
wertung gilt leider die gleiche Einschräinkung bezüglich Brauchbarkeit wie für die Cranksche 
Lôsung [12]. Weitere Verôffentlichungen von H. F. Ivey [18] sowie D. À. Bell und H. O. 
Berktay [19] brachten nur weitere Ergänzungen und aus der zusammenfassenden Darstellung 
von M. J. O. Strutt [20] ist ersichtlich, dass damit die Theorie der zylindrischen Raumladungs- 
diode noch keineswegs als abgeschlossen gelten konnte. 


2. Die Raumladungsgleichung der kreiszylindrischen Diode für konstante 
Emissionsgeschwindigkeit der Elektronen 


Angenommen seien gemäss Abb. 1 zwei konzentrisch gelegene Kreiszylinderfläichen 1 und 2 
mit den Radien rj und ro. 


CH 


Abb. 1 
Ausschnitt aus der Kreiszylinderdiode mit Emissionselektrode r1, Kontroll- 
zylinder r und Anode r2. 


Ausserdem stellen wir uns beim Radius r einen Kontrollzylinder vor. Die innere Zylinder- 
fiche 1 sei eine Kathode beliebiger Emissionsfähigkcit, der äussern Zylinderfläche 2 komme die 
Rolle der Anode zu. Den Innenraum denken wir uns evakuiert und vernachlässigen alle Neben- 
effekte wie Magnetfelder im Interelektrodenbereich, Randeffekte an den Zylindern einschliess- 
lich der Elektroden usw. Die Kathode emittiere eine gleichmässige Stromdichte j1 mit der radia- 
len Emissionsgeschwindigkeit v1 für die Elektronen. Letztere passieren dann ein Flichenelement 
des Kontrollzylinders mit der radial gerichteten, gleichfôrmig verteilten Stromdichte j mit dem 
Betrag 

2 To 1 j 1l Il ] 


(1) j= Re 
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worin 1 den gesamten Emissionsstrom bedeutet und / die Zylinderlänge, welche ohne Einschrän- 
kung gleich Eins gesetzt werden kann. 

Gefragt wird nun nach dem elektrischen Skalarpotential V des Raumladungsfeldes bei sta- 
tionärem Betrieb. Infolge der Symmetrie der Anordnung und den gemachten Voraussetzungen 
wird V eine Funktion von r allein; wir bezichen V auf den Zylinder 1 und legen das Potential 
des Zylinders 2 positiv fest: 


(2) Ft T1) —0 
AO EEE) 


Divergenz- und Energiesatz liefern folgende für den Interelektrodenraum gültige Bezichungen : 


; 164 
(3) divE y sens 
r dr 0 
(4) 2 + 
m 
à à > dd” 
worin : E = radiale elektrische Feldstärke = — 
r 
pe — Raumladungsdichte — de 
n 
L me il As 
co — Dielektrizitätskonstante des Vakuums = - - 
36710? Vm 
à dr 
1 = —— 
dé 


e — Betrag der Elektronenladung # 1,60 : 10-19 C 
m — Masse des Elektrons & 0,91 + 10-30 kg 


Unter Berücksichtigung der Beziehung (1) erhält man nach Zwischenrechnung die Gleichung 
für das Raumladungspotential 1: 


(5) 


Durch Einführen der Substitution 


€ 2/3 
(6) D] Se DEL e 
2ej1r12 m 
wird sie reduziert auf die Form: 
7 d dwW 1 
TE Ÿ = 
dr dr VW 


Der Zusammenhang zwischen der dimensionslosen Funktion W{r) und dem Raumladungs- 
potential /(r) kann auch folgendermassen angeschrieben werden: 
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D SR ANA 
WWE Eee (V + V:) 
m j1? r14 


worin Vi das der Emissionsgeschwindigkeit v1 äquivalente Potential bedeutet 


1—= He v12. 
2e 
Die Emissionsgeschwindigkeit v1 wird damit als Wirkung eines beschleunigenden, innerhalb 
des dann als vollkommen durchlässig zu denkenden Emissionszylinders 1 gelegenen Potential- 
sprunges vom Betrage 1 interpretiert. 
Je nach Anwendung wird es vorteilhaft sein, die linke Seite von (7) in einfacherer Form dar- 
zustellen. Substituiert man nämlich für die abhängige Variable r 


(8) et (:) 


/ 


so resultiert nach entsprechender Umrechnung die einfache Form für W{(x): 


dWw et 
0) AR 
dx? V W 
Diese Gleichung ist vom allgemeinen Typus 
_.. 
(10) y —— = f(x) 
F dx? 


auf welchen sich auch die Raumladungsgleichungen der parallelebenen und der kugelsymmetri- 
schen Diode reduzieren lassen. Die Funktion f(x) hat dann den Wert f(x) = 1 für den erstern 
und f(x) — (1—x) 2? für den letztern Fall. Eine interessante Untersuchung der allgemeinen 
Eigenschaften dieser Gleichung (10) gibt R. Bellman [29]; eine exakte Lôsung der Gleichung (9) 
ist jedoch bis heute nicht bekannt. Hingegen kann diese Gleichung mittels numerischer Ver- 
fahren gelôst werden; mit diesen wollen wir uns fortan befassen. 

Für den Potentialverlauf im Interelektrodenbereich der Raumladungsdiode haben wir im 
allgemeinen folgende drei Fälle zu unterscheiden : 


1. An keiner Stelle rj < r< r2 existiert ein Minimum der Potentialkurve. 
2. An einer Stelle rm tritt ein Potentialminimum mit dem Betrage Vin >> 0 auf. 
3. An einer Stelle rm stellt sich ein Potentialminimum des Betrags Vin — 0 ein. 


Bei der numerischen Behandlung der Raumladungsgleichung müssen wir uns auf die Fälle 1 
und 2 beschriänken. Die abgeleitete Raumladungsgleichung setzt nämlich die Existenz einer 
einsinnigen Elektronenstrômung für den gesamten Bereich r1 <<r<r2 voraus, während im 
Falle 3 zwischen Potentialminimum und Kathode eine doppelsinnige Strômung entsteht. Man 
zerlegt dann die Diode in 2 Raumladungsdioden rn...r1 und rm...ro je mit Emissionsgeschwindig- 
keit Null. Dies ist von anderer Seite ausführlich durchgeführt [15]. 


3. Lôsung der Raumladungsgleichung mittels Relaxationsrechnung 


Die Grundlagen, Variationen und Anwendungen des Differenzenverfahrens zur Lôsung 
partieller Differentialgleichungen sind in einigen ausgezeichneten Artikeln [21, 28] und Büchern 
[23, 25, 26] dargestellt und brauchen hier nicht erôrtert zu werden. Für die Anwendung der 
Relaxationsrechnung ist die Form (9) der Raumladungsgleichung vorzuziehen. Einmal ist dies 
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eine rein numerische Form, zum andern erhält man für konstante Schrittweite h in x gemäss 
GL. (8) eine Häufung der Punkte r gegen r1, wo im allgemeinen der Potentialgradient seine hôch- 
sten Werte erreicht, und damit einen gleichmässigeren Differenzengradienten. 

In Differenzenform kann die GI. (9) für den k-ten Punkt wie folgt angeschrieben werden: 

eTk il il 
(1) Wen + Wii —2Wg— he —— — FE Pr + — S6Wx + O (h8) = 0 
Wx 90 

Sie gilt für alle Punkte xx zwischen Null und In(r2/r1) gemäss GI. (8); 34 und 56 bedeuten die 
entsprechenden zentralen Differenzen [32]. Mit Rücksicht auf die starke Nichtlinearität der Glei- 
chung wird dann wie folgt vorgegangen: Die Randwerte W/(r1) und W{(r2) werden aus V1 und 
V2 berechnet und dazwischen eine logarithmische Potentialverteilung in r (statisches Potential 
der Zylinderdiode) angenonimen, d. h. eine lineare Verteilung für x. Für jeden Punkt wird 
dann der Wert h? Wy-1/2e?%x berechnet und daraus das entsprechende Residuum aus GI. (11) 
bestimmt. Darauf werden die Residuen mittels des einfachen (1, —2,1)-Schemas ausgeebnet, 
also wird vorderhand der Einfluss des nichtlinearen Gliedes auf die Residuen vernachlässigt. Sind 
die Residuen genüigend klein, so berechnet man neue Werte von h2W3-1/2e#%x und damit neue 
Residuen aus (11). Meist wird man mit Vorteil 34 und 36 vernachlässigen und diese Differenzen, 
wenn überhaupt, erst im weitern Verlauf der Relaxation einführen. Der jeweils entstehende Feh- 
ler kann leicht abgeschätzt werden [22, 32]. In dieser Weise ist das Relaxationsverfahren für ver- 
schiedene Randwerte von P. A. Lindsay [40] durchgeführt worden, wobei die erhaltenen Po- 
tentialverteilungen mit jenen von L. Page und N. J. Adams [17] innerhalb 8°, übereinstimmten. 

In den meisten Fällen konvergiert das Verfahren sehr rasch. Eine Ausnahme bilden alle jene 
Potentialverteilungen, bei welchen man sich im Grenzgebiet der physikalischen Gültigkeit der 
behandelten Raumladungsgleichung befindet, d. h. kurz vor Erreichen des sprunghaften Über- 
ganges zum Auftreten eines Potentialminimums mit dem Potential Null. Der Einfluss der Nicht- 
linearität der G1. (10) auf die Konvergenz des Relaxationsverfahrens ist von P. À. Lindsay [40] 
untersucht worden. Im übrigen ist jedoch über die Frage der Konvergenz bei etwas komplexeren 
und vor allem bei zweidimensionalen Problemen noch recht wenig bekannt [33, 39], auch im 
Zusammenhang mit der Netzwerkmethodel) [37]; so zieht man etwa den Schluss, dass bei Sta- 
bilität der physikalischen Anordnung der entsprechende Relaxationsprozess konvergieren müsse. 


4. Lôsung der Raumladungsgleichung mittels einer Widerstandskette 
als elektrisches Analogon 

4.1. Allgemeines 

Die Môglichkeit, mittels elektrischer Impedanznetzwerke gewühnliche und partielle Diffe- 
rentialgleichungen zu lôsen, wurde erstmals in grôsserem MaBstab in den Jahren 1939...1945 von 
G. Kron aufgezeigt. In den Zeitraum 1949...1956 fallen die zahlreichen Arbeiten von G. Lieb- 
mann über die Anwendung dieser Methode auf viele Probleme der Physik und Technik; ihm 
kommt zweifellos das Verdienst zu, diese Analogieverfahren zu dem vielseitigen und genauen 
Instrument gemacht zu haben, das sie heute darstellen. 

Mit der Lôsung der Poissonschen Differentialgleichung befassten sich zuerst D. DePackh und 
S. C. Redshaw, später G. Liebmann, R. Hechtel und G. Cremosnik [34, 36]; umfangreiche Litera- 
turangaben über diese frühern Verôffentlichungen finden sich in den zusammenfassenden Ar- 
beiten von W. W. Soroka [30] und W. J. Karplus [38]. Daneben bewies hauptsächlich G. Cre- 
moënik [34, 35, 36], dass eindimensionale, der Poissonschen Differenzialgleichung gehorchende 
Potentialfelder mittels einfacher Widerstandsketten untersucht werden kônnen. 


1) Die Konvergenz von Lôsungsverfahren bei Systemen nichtlinearer Gleichungen behandelt J. J. 
SCHÂFFER: «Contributions to the theory of electrical circuits with non-linear elements». Diss. ETH Zürich, Prom. 
Nr. 2614, (1956) und: «Zur Theorie der elektrischen Netzwerke mit nichtlinearen Elementen», Archiv f. Elek- 
trotechnik 43 (1957), 151.168. 
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4.2. Widerstandskette als Nachbildung der Raumladungsdiode 


Erteilt man der Emissionselektrode r1 das der gleichfôrmigen Emissionsgeschwindigkeit v1 
iquivalente Potential J/4, so lautet GI. (5) 


(12) dd? . 1 CF r1 j1 | 0 
dr? dr ee 2e,, 


im 


Diese Differentialgleichung muss in eine Differenzengleichung umgesetzt werden. Entsprechend 
der Feldaufteilung in Abb. 2 (oben) entwickeln wir die Potentiale Vx+1 und Vx-1 in Taylor- 


US Vk#i_ 
er, = 
vx Rk Vk Rk+1 Val 

: : Abb. 

Ik=1 lkæi Oben in der Umgebung des allgemeinen Feldpunk- 
| tes Vx mit rx = kh für die Kreiszylinderdiode. Unten der ent- 
: sprechende Ausschnitt aus der Widerstandskette mit dem Strom 
lok iox der zusätzlichen Stromquelle. 


reihen um den Mittelpunkt Vx und erhalten daraus durch Addition und Subtraktion die Ab- 
leitungen : : 


2 14 — Vry Vrx-1 — 7 D) 4 2 6 
(3) BV Vrti—Vk ne 5 Ve _d ER 
dr? h2 h2 41 dr 6! dr 
dy 14 — Vy- 1 3 de 
(14) AS A se 
dr 2h 3! dr DIN GTS 


welche wir in die GI. (12) einsetzen, um nach Berücksichtigung der Bezichung r — kh und ent- 
sprechender Zusammenfassung die Identität anzuschreiben 


V —} fl Vs —Vx il 
(15) k+1 k ( 2 | A k ee = "1 J1 E(V)=0 


2 2 
h 2k h = 11 ES 


worin F (J) den durch die hôhern Ableitungen [22, 32] bedingten Fehler darstellt: 


00 h2n d2n+1y o fjèn—-2 dny 
. ra É 
(2n + 1)! dr2n+1 “=> (2n)!  dr2n 


den wir in der Folge vernachlässigen wollen. 


Die auf die Konfiguration der Abb. 2 (unten) angewandte Kirchhofische Knotenregel liefert 
für den k-ten Maschenpunkt die Beziehung 


7 — V Vrx-1—1 
(6) Éd lo ape 
Rx+1 Rk 


die wir der mit einer noch zu bestimmenden Funktion G(r) multiplizierten G1. (15) gegenüber- 
stellen : 
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A LS PL A njihGx 


1 | 1 + 


| 2 
Tu 1 k €0 ÿ— Vx 
RE Se mienne î 
2k 2k 


Die Gin. (16) und (17) stimmen überein wenn: 


il il ; r1 1h GK 


Rk+1 ET EAST Ri = | : 10k = a 
fondée one 3 12 
— | Ci ———|G En eee 7 
2k D) k) ; 277 V 


Gx kann mit Hilfe der Eindeutigkeitsbedingung für die Widerstände bestimmt werden. So muss 
Rx als «unterer» Widerstand der Stelle (k+1) identisch sein mit Rx+1 als «oberem» Widerstand 
der Stelle k: 


1 Î 
1 ÿ 1 
ebrihe lle 
| mn) k+1 | | k 
d. h. 
Gk+1 _ Gk Gi 
k+1 k 


Mit Gr = kG1 und k — 1 resultiert aus der Formel für Rx der Widerstand R1 — 2/G1, den wir 
als Bezugswert wählen; die Reïhe lautet alsdann : 
R: R; Ri 


(18) Ri , TRS nn MAN O0 LD , Rx le — — 
5 5 2k—1 


während der Strom iox den Betrag 


le) . | 
(19) _— AU €] 


Te = —— 

2e VX 

£o Ri || _ Vr V F 
NI 


aufweist. Damit sind alle für die Feldnachbildung notwendigen Grüssen festgelegt. 

Für die Bestimmung eines Raumladungsfeldes bei gegebenen Abmessungen r1 und ro, 
Emissionsstromdichte j1, Emissionsgeschwindigkeit v1 — V11/2 und Anodenpotential V2 muss 
nun genau wie bei der Relaxationsrechnung iterativ vorgegangen werden. Zur Berechnung der 
iox als Nachbildung der Raumladung sollte nämlich neben j1 das Potential V4 bereits bekannt 
sein. Die Stromdichte kann auf den Langmuirschen Wert bezogen werden, während wir für das 
Potential das bekannte statische Feld x — V2 : In (rx/r1) / In (ra/r1) für pe — 0 als nullte Nähe- 
rung benutzen (d. h. die Lôsung der Laplaceschen Gleichung AV — 0). Damit lassen sich die iox 
einstellen und mit der resultierenden 1. Näherung des Potentialverlaufes dann neue 0x berechnen. 
Dies wird so lange fortgeführt, bis sich die V4 nicht mehr ändern. 

Die Diodenstrôme betragen nach I. Langmuir für innenliegende Kathode bei der Emissions- 
geschwindigkeit v1 — 0 


Vo3/2 
(20) Hoi") 
r1 ro P?10 


bzw. für aussenliegende Kathode 
V,3/2 
(21) ln == 


Or) 
r2 11 4221 
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wWorin 


a — Diodenkonstante — (49/9) (2 e/m)1/2 — 2,34 + 1076 AV-3/2 


8219 — Langmuirsche Funktion 8? für innenliegende Kathode?), 6219 — 82 : ; 


Beide Funktionen sind in Abb. 9 für den Bereich 1  r2/r1 5 aufgetragen. 


4.3. Messungen und Auswertung 


Für die Messung der Raumladungskennlinien wurde die Anordnung nach Abb. 3 benutzt. 
Der Bezugswiderstand für die Kette beträgt R1 — 4,5 kQ, alle Kettenwiderstände sind auf 2%, 
genau abgeglichen. Die für die gegenseitige Nichtbecinflussung bei der Einstellung der iox not- 
wendige Stromquelleneigenschaft wird durch Reihenschaltung von 5,5 MQ mit der stabilisierten 
Speisespannung 500 V gewährleistet, wobei 5,5 MQ — 5 MQ (Grob-) + 0,5 MQ (Feinregelung) 
ein gewühnliches Kohle-Doppelpotentiometer ist. Statt der Kombination Ampèremeter + Ein- 
stellpotentiometer kôünnten auch geeichte Einstellwiderstinde verwendet werden [40], doch ist 
das notwendige Auflôsungsvermügen nur mit Wendelpotentiometern zu erreichen; eine gegen- 
über der gewählten sehr teure Lüsung. Interessant ist die automatische Einstellung der iox mittels 
Analogcomputer-Einheiten [31], doch wird der grosse Aufwand pro Zweig (3 Rechenverstärker, 
1 Servoverstärker, 1 Motor, 3 gekuppelte Wendelpotentiometer und entsprechende Speisung) 
diese Lôüsung wohl selten gestatten. 


500 V 
A 
VA V, 
0.10 V 
Abb. 3 Ausschnitt aus der zur Messung verwendeten Wider- 
standskette mit grundsätzlicher Messanordnung. Die 
} Voltmeter mit dem Innenwiderstand R; messen Ua — 
Vo — Vi und Ux = Vrx — V1. Es muss Ri > Rx. 


Die Kathode erhält das Potentional V1, die Anode V2. Ein Voltmeter kontrolliert die K-A- 
Spannung Ua — Vo — V1, welche mit der Speisespannung 500 V vom gleichen Hochkonstant- 
Netzgerät (Sfeinlein HK 101) geliefert wird. Gemessen wird das Potential Vx = V1 + Uk. 

Nachgebildet wurde eine Diode mit r2/r, — 5 und innenliegender Kathode, wobei r1 = 4h 
und r2 — 20h; dies ergibt 15 iok-Zweige zur Nachbildung der Raumladung. Für die Auswertung 
ist es von Vorteil, die Emissionsstromdichte j1 auf j12 nach Gl. (20) zu beziehen : 


2) NN Be. 


j12 a V23/2 


2) Der erste Index bedeutet die Kathode, der zweite die Anode, wobei immer r9 > rm > 1. 
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so dass sich aus Gl. (19) die ioK berechnen zu: 


. 2 a V93/2 C19:K 
(23) IQ —— a 


He 7 
co R1 B212 : e Van V k 


Für die gewählte Diode ergibt dies bei V2 — 10 V 


Vx 


Die Raumladungskennlinien wurden nun für den Bereich 0 < (V/1/V2) 1 in Schritten von 
0,1 für je etwa zehn «-Werte bestimmt; eine typische Kurvenschar zeigt Abb. 4 für V1/Vo = 0,7. 
Daraus wurde unter Zuhilfenahme hôherer Differenzen und Interpolation [27] der Ort rm/r2 
und der Betrag Vin/Vo des Potentialminimums bestimmt. In kritischen Fällen kann auch die 


Abb. 4 

Typischer Verlauf der Potentialkurven (1//V2) = f(r/r1), 
hier für V1/V2 = 0,7. Die oberste Kurve gibt das statische 
Feld ohne Raumladung, die folgenden 10 Kurven das 
Raumladungspotential für die «-Werte 0,75...1,00...1, 
50...2,00...2,50...3,00...3,50...3,75...4,00...4,20. Die einzel- 
nen Potentialminimumspunkte sind verbunden. P2 ist ein 
Punkt der Grenzkurve 2. 


Abb.5 Kurvenschar für (V/m/V2) = f(K) mit 
(V/V2) als Parameter in Schritten 
von 0,1. Das «Minimumsgebiet» 
wird durch die Grenzkurven 1 und 
2 eingeschlossen; die Grenzkurve 3 
fällt mit der «k-Achse zusammen. 
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Lüsung mittels der Kette als nullte Verteilung des Relaxationsprozesses betrachtet werden und 
mittels des im vorigen Abschnitt behandelten Verfahrens verbessert werden. 


Die bezogenen Grüssen Vm/V2 und rm/re sind in den Abb. 5 und 6 in Funktion von K mit 
Vi/Vo als Parameter aufgetragen, während Abb. 7 den Zusammenhang zwischen V1/V2 und 
Vm/V2 mit k als Parameter darstellt. 


25 


0 | 2 3 4 5 — 


Abb. 6 Kurvenschar für (rm/r1) — Ê (K) bzwW. (rm/r2) = f (x). Die Kurven für V1/V2 = 0,2 und V1/V2 — 0,1 sind der Deutlich- 
keit halber nicht eingetragen, da hier infolge des Verlaufs der Grenzkurve 1 (s. Abb. 5) Überschneidungen auftreten. 
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—, 
Abb. 7 Kurvenschar für (Vm/V2) = f (V1/V2) mit der normalisierten Stromdichte als Parameter. 
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4.4. Grenzkurven 
Besondere Aufmerksamkeit beanspruchen folgende 3 Grenzkurven: 


GK 1: Geometrischer Ort der kleinsten Werte r (11/12), für welche im Interelektrodenraum 
ein Spannungsminimum Vm >> 0 mit rm >> 1 existiert, d. h. GK zwischen den Fällen: 
PR tn TU OV TT Tu 1. 

GK 2: Geometrischer Ort der grôssten Werte k (V1/V2), für welche im Interelektrodenraum 
ein Spannungsminimum besteht, d.h. sprunghafter Übergang vom Zustand 0<V<11, 
fl fm < te zum Zustand Vin = 0, f] < tm < re. 


GK 3: Geometrischer Ort der Werte k (11/12), für welche sich im Interelektrodenraum gerade 
ein Potentialminimum mit Vn = 0, #1 < rm < r2 ausbildet. 


Der Übergang von der Grenzkurve 1 zur Grenzkurve 2 ist kontinuierlich und umfasst das 
gesamte «Minimumsgebiet» (ausgezogen in Abb. 5), jener von Grenzkurve 2 zur Grenzkurve 3 
ist sprunghaft (gestrichelt in Abb. 5). 

Die Bestimmung der Grenzkurve 1 aus den gemessenen Raumladungskurven ist ohne 
weitere Verkleinerung der Maschenweite nicht sehr genau. Es wurde deshalb zu folgender 
Methode gegrifflen: Für r’1 < r1 erhalten wir eine Diode r’1...r2, für welche bei gegebenem 
V1/Ve und einem bestimmten Wert « gerade r'm — r1 ist; das Wertepaar (k, V1/V2) ergibt dann 
einen Punkt der Grenzkurve. 

Die Festlegung der Grenzkurve 2 ist einigermassen mühsam, da das Verfahren hier schlecht 
konvergiert, so dass die Anzahl der Iterationsschnitte (& 4 im Minimumsgebiet) hier verdoppelt 
bis verdreifacht werden muss. Einige Punkte wurden deshalb mittels des — auch nicht besser 
konvergierenden! — Relaxationsverfahrens kontrolliert. 


f Fm Lg) f M Tm r 


Abb. 8 Zur Berechnung der Grenzkurve 3. Zerlegung der Diode r1—r2 in die Dioden rm—r1 und rm—r2 bei Vorliegen eines 
Potentialminimums Vm = 0 für r = rm. 


Demgegenüber kann die Grenzkurve 3 in einfacher Weise berechnet werden. Nach Abb. 8 
wird die Diode r1...ro zerlegt in die zwei Dioden rm...r1 (Raumladung Z Strom 2 11 — 12) und 
Fm….r2 (Raumladung © Strom 12); für beide gilt dann das Langmuirsche Gesetz: 


V13/2 
(24) 2 1 — Le — ne (27 fini) 
Fm 1 X?m1 
| V53/2 
5) = <— 27 rm l) 


Aus diesen zwei Beziehungen kann das gesamte Gebiet, in welchem ein Potentialminimum 
Vm = 0 besteht, berechnet werden [15]. 


ol 


Setzen wir nun 2 11 — 19 — Jo und normieren mit 112 nach GI. (20), so bestimmen die zwei 


Gleichungen 
; 0,7666 
(26) B?m2 D 
K 
V K à 2/3 
(27) — É | 
Vo 9599 


den Verlauf der Grenzkurve 3; der Zusammenhang zwischen B2m2 und «2h ist dabei der Abb. 9 
zu entnehmen. Es ergibt sich, dass für V1/V2 > 0,1 die Grenzkurve mit genügender Genauigkeit 


durch die Gerade 
(28) 


"i 


Kk= 10614 |— 


Vo 


approximiert werden kann. Die Grenzkurve 3 ist zusammen mit den Grenzkurven 1 und 2 als 


k (V1/V2) in Abb. 10 dargestellt worden. 


m2 mi 
10 10 
09 9 
08 8 
07 7 
0,64 L6 
0,5 | 5 
04 | 4 
031 3 
02 2 
01 1 
0 0 
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Abb.9 Langmuirsche Korrekturfunktionen (22 (rAnode/"Kathode) bei innen- 
liegender Kathode und «2 (rxatnode/rAnoae) bei aussenliegender 
Kathode für den Bereich 1 < rofrm < 5 für B2m2 bzW. 1 < rm/ 
r1 < 5 für œny nach [9]. 
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Abb. 10 Verlauf der Grenzkurven 1,2 und 
3 als k (V1/V2). 
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